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Abstract 



Results are reported from a search for the anomalous production of highly boosted Z 
bosons with large transverse momentum and decaying to the dimuon final state. Such 
Z bosons may be produced in the decays of new heavy particles. The search uses pp 
collision data at a/s = 7 TeV, corresponding to an integrated luminosity of 5.0 fb^^ 
recorded with the CMS detector. The shape of the observed transverse-momentum 
distribution of Z bosons is consistent with standard-model expectations. Constraints 
are obtained on models predicting the production of excited quarks decaying via elec- 
troweak processes. For excited-quark decays involving only standard-model gauge 
bosons and coupling strengths, masses of excited quarks below 1.94 TeV are excluded 
at the 95% confidence level, assuming a compositeness scale equal to the excited- 
quark mass. For excited-quark production via novel contact interactions, masses be- 
low 2.22 TeV are excluded, even if the excited quarks do not couple to gluons. 
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In the standard model (SM) of particle physics, the transverse-momentum (px) spectrum of Z 
bosons produced in high-energy pp collisions is predicted to be a smoothly falling distribution. 
A broad range of new physics models such as quark compositeness [1 . 2J, supersymmetry [3J, 
technicolor [4J, and extensions of the standard model with new gauge groups predict de- 
cays of heavy particles involving Z bosons, which would introduce deviations from a smooth 
distribution of the Z transverse-momentum spectrum. The Z boson transverse-momentum dis- 
tribution was measured by ATLAS [6| and CMS [7J in pp collisions at i/s = 7TeV using about 
40 pb^^ of integrated luminosity, showing good agreement with next-to-next-to leading order 
perturbative QCD calculations. This Letter presents a study of the inclusive dimuon transverse- 
momentum spectrum, pj{mi), of Z boson production in pp collisions. The measured spectrum 
is used to search for a new heavy particle decaying into high momentum Z bosons. Because the 
measurement is inclusive, no constraints are imposed on the presence of additional particles 
in the decay of the hypothetical heavy particle. The integrated luminosity is also significantly 
larger than that used in previous studies [jH |Zl- The final results are interpreted within the 
framework of specific models of excited-quark production, q* qZ, Z ^^H^ ■ 

Compositeness models explain the observed mass hierarchy of quarks and leptons by intro- 
ducing quark constituents, predicting a multitude of excited-fermion states, including excited 
quarks mill. In the effective Lagrangian describing the gauge-mediated transitions of excited 
fermions, the couplings to the strong and electroweak sectors are measured in units of the 
strengths of the SM gauge couplings gs (strong coupling), g = e sin 6yi, and g' = e cos 6w, where 
e is the electron charge and Oyj is the weak mixing angle. Thus, the corresponding strengths 
for the new interactions are parametrized in terms of scale factors fs,f, and /' with respect to 
the SM couplings. In many studies, new couplings are assumed to be equal to the SM gauge 
couplings (i.e. f = f = fg = 1) and the compositeness scale. A, is usually taken to be equal 
to the excited-quark mass, Mq* . Moreover, the excited-quark model gauge theory can be ex- 
tended with novel four-fermion contact-interaction terms arising from new strong d5mamics. 
Excited quarks are usually sought in the dijet final state [8J in hadron collisions, but with /s = 
there is no sensitivity. Gauge- and contact-interaction transitions are also typically assumed for 
excited-fermion searches at HERA |9 1, while contact-interaction production is generally probed 
in excited-lepton searches at hadron colliders [10|. 

Searches for q* production use a range of different strategies. The DO experiment searched for a 
mass resonance in the Z plus jet system, with the Z boson detected via its dielectron decay mode 
HU. The search presented in this Letter tests the spectrum for signatures of boosted Z 

boson decays inclusively, without specifying the recoiling system so as to be less model depen- 
dent. A reduced parameter space of the excited-quark production models is probed assuming 
/ = /' = ! and using three independent parameters: the mass of the excited quark, Mq*, the 
compositeness scale. A, and the strong coupling scale factor, fg. Only first generation excited 
quarks, degenerate in mass (u*, d*), are considered, as they have the largest production cross 
sections in proton collisions. The u* (d*) branching fraction to Z is 3% (5%) for gauge-only 
couplings with fg = 1, increasing to slightly over 20% (30%) with fg = JlUZl. The production 
of excited quarks via a gauge transition (qg — > q*) is treated separately from the production 
via a contact-interaction diagram (qq' — > q"q*), though in both cases only the gauge decay 
to a Z boson (q* — > qZ) is considered. When contact-interaction production is assumed, such 
terms are also included into the Z branching fraction calculation. The gauge-interaction tran- 
sitions imply that Z bosons arise from two-body final states (qZ), whereas contact-interaction 
production yields Z bosons in three-body final states (qqZ), as illustrated in Fig.[l] 

This study uses proton-proton collisions data a\. ^Js = 7 TeV, corresponding to an integrated 
luminosity of 5.0 fb^^ recorded by the Compact Muon Solenoid (CMS) experiment at the Large 
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Figure 1: Feynman diagrams for the production of q* —> qZ, Z — > via (left) a gauge 

interaction and (right) a contact interaction. Charge conjugation of an initial state quark is 
implied in both diagrams, as is permutation of the quark isospin in (right). 

Hadron Collider (LHC). The central feature of the CMS apparatus is a superconducting solenoid, 
of 6 m internal diameter, providing a field of 3.8 T. Inside the magnet coil are the silicon pixel 
and strip tracker, the lead tungstate crystal electromagnetic calorimeter, and the brass/ scintillator 
hadron calorimeter. Muons are detected in gas-ionization detectors embedded in the magnet 
steel return yoke. In addition to the barrel and endcap detectors, CMS has extensive forward 
calorimetry. The trigger system, composed of a custom hardware layer feeding into a commer- 
cial processor farm (High Level Trigger, HLT), reduces the event rate to approximately 300 Hz 
for storage and further analysis. A detailed description of the CMS apparatus may be found 
elsewhere |T2|. 

Anomalous production of Z bosons arising from heavy new particles and decaying into dimuon 
final states is characterized by a pair of oppositely charged isolated muons with an invariant 
mass consistent with that of the Z boson and a high dimuon transverse momentum. The dom- 
inant irreducible background is due to SM Z — ?► ji^fi^ production. The other background 
sources considered in the analysis are from prompt muon processes. Prompt muons are de- 
fined to be muons originating directly at the primary vertex or originating from the decays 
of short-lived particles. The background sources of prompt muons considered are tt, diboson 
(WW, WZ, ZZ), and Z — > t+t^ — )■ X production. In addition, jets may be misidenti- 

fied as muons and contribute to the dimuon transverse-momentum spectrum through multijet 
and W plus jets final states (non-prompt muon backgrounds). The backgrounds were mod- 
eled using simulated samples produced with the full Geant4 [13] based CMS detector sim- 
ulation, trigger emulation, and event reconstruction chain. Different samples of SM Z boson 
production were generated with POWHEG vl.l I[l4l - [T6) and the MadGraph matrix element 
generator ||T7[. both interfaced to the PYTHIA v6.424 IflSl parton shower generator. Diboson 
and QCD processes with a muon in the final state were modeled using pythia. Events from 
tt and W plus jets were modeled using MadGraph interfaced to the pythia parton shower 
generator. For the excited-quark modeling, we rely on simulations assuming separately either 
gauge- or contact-interaction production, generated for mass points ranging from Mq* = 500 
to Mq* = 2300 GeV in steps of 100 GeV using pythia. Gauge decays to a Z boson are assumed 
for both production choices. 

Events are selected offline to have two high-px (Pt > 45 GeV), oppositely charged, isolated 
muons. Events used in the analysis were collected using a single-muon trigger. The algorithm 
requires a muon candidate to be found in the muon detectors by the first level trigger system. 
The candidate track is then matched to a silicon tracker track, forming a HLT muon. The HLT 
muon is required to exceed a pj threshold of 40 GeV and to be reconstructed in | < 2.1, where 
the pseudorapidity is defined as // = — ln[tan(0/2)], where 9 is the polar angle with respect to 
the direction of the counterclockwise beam. A right-handed coordinate system is used in CMS, 
with the origin at the nominal collision point, the x axis pointing to the center of the LHC ring, 
the y axis pointing up (perpendicular to the LHC plane), and the z axis along the anticlockwise- 
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beam direction. The azimuthal angle (p is the angle relative to the positive x axis measured in 
the x-y plane. 

Muon candidates are reconstructed with two algorithms fl^. In the first "tracker-muon" algo- 
rithm, tracks are fit with hits in the silicon tracker, propagated outward, and matched to hits 
in the muon system. In the second "global muon" algorithm, a global fit is performed to hits 
both in the silicon tracker and the muon system. At least one muon candidate is required to 
be successfully reconstructed by both algorithms. Requiring both muons to be reconstructed 
as global muons would eliminate some background, but would introduce an efficiency loss 
associated with muon pairs with small opening angles, AR = \/ A(p'^ + Af]^ < 0.3, which is 
typical of highly boosted Z candidates. The track associated with each muon candidate is re- 
quired to have hits in at least 8 layers of the silicon tracker, at least one hit in the pixel detector, 
and a magnitude of the transverse impact parameter below 0.2 cm to be consistent with a par- 
ticle emanating from the primary interaction vertex. The candidate reconstructed as a global 
muon is further required to have hits in at least two different muon detector stations llT9l . The 
momentum vector of the offline muon candidate reconstructed with the global fit algorithm 
must be matched in direction to the HLT candidate that triggered the event. If both muons 
are reconstructed with the global-fit algorithm, only one of them is required to match an HLT 
muon satisfying the trigger requirements. Both muons must have pj > 45GeV and \ f]\ < 2.1. 
The trigger efficiency for a single muon that passes all offline selection criteria is 92%. Muon 
candidates from the same vertex are selected by performing a common-vertex fit and requiring 
the vertex to be below 10. The muon pair is required to have an invariant mass consistent 
with the Z boson, 60 GeV < M,,,, < 120 GeV. 

Finally, events are required to contain at least one reconstructed primary interaction vertex with 
at least four tracks, located within 2 cm of the center of the detector in the direction transverse 
to the beam and within 24 cm in the direction along the beam. These requirements help to 
reject cosmic-ray background. Additional suppression of cosmic-ray muons is obtained by 
requiring the three-dimensional opening angle between the two muons to be smaller than n — 
0.02 radians. 

For the production processes considered here, simulations show that the muons from Z decay 
are usually isolated from the hadronic activity in the event. To measure the isolation for each 
muon candidate, a cone is constructed of radius AR < 0.3 around the track direction at the 
primary interaction vertex. We estimate the activity around the lepton by computing the scalar 
sum of the transverse momenta of charged tracks within this cone, excluding the muon candi- 
date itself. Only tracks with distance-of-closest approach to the primary vertex less than 0.2 cm 
along the beam direction are considered. The ratio of the this isolation quantity to the px of 
the muon itself is required to be less than 0.1. The efficiency of this requirement for dimuon 
events in the Z mass window is found to be independent of the number of reconstructed pri- 
mary vertices, showing that the isolation is insensitive to the effects of multiple pp interactions 
(pileup), which is on average nine per bunch crossing. The simulated samples are re-weighted 
such that the distributions of the number of reconstructed primary vertices correspond to the 
measured distributions in data. For highly boosted Z boson decays, the two muons can ap- 
proach collinearity, and therefore can affect each other's isolation calculation. We correct the 
isolation measurement of a muon by excluding the other muon candidate from the isolation 
sum if the two are closer than AR = 0.3 and if the scalar sum of transverse momenta defining 
the isolation before the correction is applied is greater than 90% of the transverse momentum 
of the other muon. 

The distribution of the inverse of the observed dimuon transverse momentum (1/ pj^fiji)) is 
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shown in Fig. |2]for candidates passing all the selection criteria and having 1/ pj{fiji) less than 
0.008 GeV"^ (Pt(Fj") > 125 GeV). The 1/pt(/^jM) distribution compresses the salient features of 
the pj[]dji) distribution onto an abscissa of finite range. There are 7044 (29) events observed 
for 1/ fj^iji) less than 0.008 (0.002) GeV^^. The dimuon candidate with the highest transverse 
momentum is found to have Px(mf) = 940 GeV. 

All background components other than Z — > ji^ji^ production contribute less than 2% and 
have negligible impact on the overall shape of the 1/ pi{}i}i) spectrum. The simulation pre- 
dictions were cross-checked by estimating the prompt and non-prompt muon background 
contributions from the data. The dominant non-Drell-Yan background contributions at high 
dimuon transverse momentum arise from tt and diboson production. All these processes are 
flavor symmetric and produce twice as many e}i pairs as ee or jiji pairs. The prompt muon 
background from tt is estimated by comparing the transverse-momentum spectrum of e^}V^ 
events between data and simulation and correcting for differences in the geometric accep- 
tances and efficiencies. This method predicts 46 (5.3) isolated pairs with \/ipj{}i}i) < 
0.008 (0.005) GeV^^. The estimate of the residual contribution from background events with at 
least one non-prompt or misidentified muon was made by looking at events selected from the 
data sample with single muons that pass all selection cuts except the isolation requirement. A 
map is created with probability of a muon to pass the isolation criteria employed in the anal- 
ysis as a function of px arid r]. This probability map is corrected for the expected contribution 
from events with single prompt muons from tt and W decays and for the observed correlation 
between the probabilities for two muons in the same event. The probability map is used to pre- 
dict the number of background events with two isolated muons based on the sample of events 
that have two non-isolated muons. This procedure, which has been validated using simulated 
events, predicts a background of less than one event for 1/ pj(jifi) < 0.008 GeV^^. More details 
on similar techniques are described in Ref. ||20| . 

The main background contribution, which arises from SM processes producing Z — > , 
is evaluated from a template fit to the observed data distribution for events with 1/ pj{}i}i) < 
0.008 GeV^^. The choice of the analytical template is driven by studies with Z — > simu- 
lated events generated with POWHEG and MadGraph and hadronised using PYTHIA: 

y(;c) =erf(fl-x^-c)-erf(-c), (1) 

where erf is the error function and x is 1/ pj{}i}i). The parameters a, b, and c of the background 

analytical template are obtained by fitting the region 1/px(^'M) ^ [0.0028,0.008] GeV^^. This 
choice avoids potential contamination from a signal with Mq* > 1.0 TeV, which could oth- 
erwise bias the fit. This effect has been tested with Monte Carlo pseudo-experiments based 
on analytical forms of the 1/ pj{j.iji) spectrum to describe both the background and the signal 
shapes. Conversely, the number of events in the lower 1/ px(f<H) region provides a robust esti- 
mation of the background, hence potential contamination due to q* signals with Mq* < 1.0 TeV 
negligibly biases the fit. Simulated Drell-Yan events were used to demonstrate that the number 
of events predicted by the fit to the 1/ pj{}i}i) distribution is not sensitive to the choice of the 
1/ pj{nji) fitting region. The resulting fit to the data is shown in Fig.|2] 

The observed 1/ pj{}ifi) spectrum shape agrees with expectations based on SM processes. With 
no clear evidence for new physics, we proceed to set 95% confidence level (CL) upper limits 
on the cross section (a) for an excited-quark production and decay process q* — > qZ. For 
the calculation of the limits, we adopt the frequentist construction CLg ||2T1 122| . We use an 
extended likelihood Il23| , built from characteristic signal and background probability density 
function templates. In all cases, we use the RooStats implementation of the algorithm [24]. 
In order to parametrize the underlying shape from the excited-quark decay model, we rely 
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CMS, 5.0 fb"^ at ^=7 TeV 




1/p>^i)[GeV-^] 

Figure 2: Distributions of \/pj{}i}i) for data (points with error bars), simulated SM back- 
grounds (stacked histograms), and simulated signal models (overlaid histograms). The signal 
models use either gauge interactions (GI) or contact interactions (CI). The blue solid line corre- 
sponds to the analytical template fit in the region defined as l/pT{}iji) G [0.0028, 0.008] GeV^^. 
The total number of SM background events is rescaled to the number of events observed in 
data. The normalization of the signal distributions is arbitrary. 

on simulations of models assuming only gauge-interaction production and models assum- 
ing only contact-interaction production. The upper limit on the cross section for a model 
of interest is obtained from the upper limit on the number of signal events divided by the 
integrated luminosity, J Cdt, and the detector acceptance times the efficiency. Axe. The 
acceptance. A, is defined as the fraction of generated dimuon candidates that have invari- 
ant mass 60GeV < M,,,, < 120 GeV and both muons with \?]\ <2.1 and pj > 45GeV. For 
the excited-quark model produced with gauge interaction, the signal acceptance times effi- 
ciency after the complete selection criteria for dimuon l/pj{jiji) < 0.008 GeV^^ is between 
42% (Mq* = 0.5 TeV) and 73% (Mq* = 2.0 TeV). For the contact-interaction production scenario, 
the acceptance times efficiency is between 51% (Mq* = 0.5 TeV) and 71% (Mq* = 2.0 TeV). 

The values for the Axe are obtained from simulation, but the efficiency of passing the offline 
selection criteria, e, is corrected by a factor Se, which accounts for data-simulation discrepan- 
cies. The single-muon triggering and particle identification efficiencies are measured in a sam- 
ple of inclusive Z — >• events in data and MC simulation separately using a tag-and-probe 
technique IIT9ll25ll . The difference between the efficiency measured in data and simulation is 
used to correct the dimuon selection efficiency, e. The final value of Se used in the analysis is 
0.98 ± 0.03. Most of the systematic uncertainty in the scale factor arises from the low-statistics 
region where AR between the muons is <0.3. We assign a systematic uncertainty of 2% for 
the detector acceptance, evaluated by variations of final- and initial-state radiation, as well as 
parton distribution functions (PDF) sets in simulations of the signal models. The PDF uncer- 
tainties on the final selection signal acceptance have been calculated using the PDF4LHC 126] 
prescription, using the PDF sets CT [27], MSTW |28| and NNPDF |29|. The PDF uncertainties 
are found to be well below 1% for the q* model considered. The systematic uncertainty in the 
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integrated luminosity is estimated to be 2.2% ||30| . The background yield is treated as a nui- 
sance parameter in the fit procedure. Its systematic uncertainty comes directly from the data 
fit and is estimated to be 2%. The uncertainties derived on the shape parameters of the tem- 
plate fit do not affect the final results. Other fit templates were tried, but gave fits that tended 
to over-estimate the background yield, resulting in more stringent limits. For each systematic 
uncertainty, we use a log-normal distribution for the nuisance parameters in the likelihood 
construction. 

The 95% CL upper limits on cross section for both gauge interaction and contact interaction 
are shown in Fig.|3] The expected limits and the one- and two-standard-deviation limit bands 
are overlaid on top of the observed limit. Generally, the limits are within the two standard 
deviation limit bands. A small deviation is found in the region starting from \/pj{}i}i) ~ 
0.002 GeV^^ resulting in higher limits for the gauge-interaction models with Mq* G [1.0-1.4] TeV. 
Since the gauge interaction yields a two-body final state while the contact-interaction produc- 
tion yields a three-body final state, the gauge-interaction signal distributions are generally nar- 
rower than the contact-interaction distributions. Owing to this and the sharp left edge of the 
gauge-interaction distributions (Fig. |2]|, a small excess is favored in this region. The limits at the 
q* mass value probed are correlated as a result of the large width of the signal shapes present 
in the l/pj{}i}i) spectrum. 
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Figure 3: 95% CL upper limits on the cross section of q* — > qZ for the benchmark models 
derived with the CLg limit setting criterion: (left) gauge-interaction production, (right) contact- 
interaction production. The results shown are obtained assuming all new couplings equal to 
the SM couplings between ordinary fermions (/ = /'= /s = 1) and A = Mq*. 

The limits on cross section are then translated into contours in the parameter space of the 
gauge-interaction and contact-interaction models. For a given pair of values of fs and Mq*, 
we probe different values of the compositeness scale A in steps of 100 GeV and select the value 
of A producing the closest cross section to the observed limit. We use pythia to compute 
the u* and d* cross section values at leading order. The excited-quark branching fractions are 
computed according to the equations reported in Ref. HI. In Fig.|4]we report the final limit con- 
tours in the Mq* and A plane for different assumptions of the strong coupling, fs , separately 
in gauge-interaction and contact-interaction production. The limit contours directly translate 
the results obtained on the cross section: a lower value for the upper limit on the cross section 
corresponds to a higher value for the lower limit on A; e.g., the downward statistical fluctua- 
tion at Mq* = 700 GeV in the gauge-interaction production model translates into a bump in the 
corresponding A contour. 

In summary, a search for anomalous production of highly boosted Z bosons in the dimuon 
decay channel has been performed using proton-proton collision data with an integrated lumi- 
nosity of 5.0 fb^^ collected at ^^s = 7 TeV by the CMS experiment. The Z transverse-momentum 
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CMS, 5.0 fb"^ at\fs=7 TeV 
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Figure 4: Contours of excluded parameters for the gauge- and contact-interaction q* produc- 
tion mechanisms. For a given fs assumption, each (Mq*, A) below the limit curve is excluded. 
The bump in the A contour at Mq* = 700 GeV for the gauge-interaction production models 
is a direct consequence of the downward fluctuation on the upper limit on the cross section 
observed at Mq* = 700 GeV (Figure |3|. 



distribution observed is consistent with SM expectations. Limits are derived on the specific 
model of excited-quark production and decay in q* — > qZ. We report 95% exclusion contours 
in the compositeness scale versus excited-quark mass plane for two production scenarios and 
for several choices of the relative coupling to gluons. Under the assumptions for the param- 
eters Mq* = A and / = /' = /s = 1/ our limits exclude excited quarks at 95% CL with 
Mq* < 1.94 TeV for gauge production and Mq* < 2.15 TeV for the contact interaction respec- 
tively. In comparison, the corresponding exclusion for the gauge production of q* from a search 
in the dijet final state is Mq* < 2.49 TeV IH. For q* created through the contact-interaction sce- 
nario, we exclude a large section of the strong coupling phase space, including /s = with a 
mass up to 2.22 TeV, much higher than the /s = limit of 252 GeV set by HI [9J in ep collisions. 
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